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ABSTRACT 



We present Submillimeter Array observations of the young stellar object IRAS 04579+4703 in the 1.3 mm continuum and in the 
I2 CO(2 - 1), 13 CO(2 - 1) and C 18 0(2 - 1) lines. The 1.3 mm continuum image reveals a flattened structure with a mass of 13 M Q . 
The 12 CO(2 - 1) line map and position-velocity (PV) diagram, together with the broad wing (full width = 30 km s~') of I2 CO(2 - 1) 
line, clearly show that there is an outflow motion, which originates from an embedded massive YSO in this region. The lengths of 
the blue-shifted and red-shifted lobes are 0.14 pc and 0.13 pc respectively. The total gas mass, average dynamical timescale and 
mass entrainment rate of the outflow are 1.8 M Q , 1.7 x 10 4 yr and 1.1 x 10~ 4 M Q yr , respectively. The flattened morphology of 
the continuum source perpendicular to the outflow direction, and the velocity gradient seen in the spectra of C 18 0(2 - 1) taken from 
different locations along the major axis of the continuum source, suggest the presence of an accretion disk in this region. 

Key words. ISM: individual (IRAS 04579+4703) — ISM: kinematics and dynamics — ISM: molecules — stars: formation 



1. Introduction 

Observational evidences suggest that massive stars may be 
formed through the accretion-disk-outflow processes similar to 
low-mass stars (e.g, Shu et al. I1987I Cesaroni et al. 1999; 
Shepherd et al. 2001;Patel et al. [20051 Jiang et al. [20031 
Zinnecker & Yorke 2007). Searching for collapse, accretion and 
outflow with interferometers in massive star formation regions 
has boomed (e.g, Beuther e al. I2004t Qin et al.[2008; Furuya et 
al. 1201 11 1, but the limited observations cannot provide detailed 
information for our understanding of massive star formation. 
The major difficulty is due to the large distances (> 1 kpc), clus- 
tered formation environments and shorter evolutionary timescale 
of massive stars. 

IRAS 04579+4703 is a young high mass stellar object 
(YSO). Molinari et al. d 19961 1 carried out the NH 3 observa- 
tion toward IRAS 04579+4703, and obtained a distance of 2.47 
kpc from the sun and a bolomenic luminosity of 3.91xl0 3 L Q . 
However, Molinari et al. (1998a) failed to observe any 6 cm 
emission from IRAS 04579+4703. Sanchez-Monge et al. ([20081 
performed higher spatial resolution continuum observations at 
1.2 mm, 7 mm, 1.3 cm and 3.6 cm toward IRAS 04579+4703. 
Their observations showed that the 7 mm emission, similar to 
the 1.3 cm and 3.6 cm emission in morphology, is elongated in 
the southwest-northeast direction. The location of their 1 .2 mm 
continuum peak agrees well with that of the source detected at 
cm wavelengths as well as an H2O maser (Palla et al. 119911 & 
Migenes et al. 1999), indicating that the emissions only concen- 
trate on a small and compact region and the IRAS 04579+4703 
is at an early evolutionary stage. Using the SED fit, Sanchez- 
Monge et al. ( 120081 derived a dust temperature of about 30 K 
and obtained a gas mass of 23 M G . The line 12 CO(2 - 1) obser- 



vations of Zhang et al. (2005) did not reveal molecular outflow 
from IRAS 04579+4703, while Wouterloot & Brand ( fT989l ) ob- 
served the 12 CO(l - 0) emission with a linewidth of 6. 1 km s . 
Additionally, Varricatt et al. (2010) detected jet-like H2 knots in 
the northwest-southeast direction. However, these observations 
can not provide sufficient evidences for the outflow. Higher spa- 
tial resolution observations of the molecular lines and continuum 
in submillimeter/millimeter waveband are needed to reveal the 
detailed source structure and kinematics in this region. 

In this paper we present Submillmeter Array (SMA) obser- 
vations of the young stellar object IRAS 04579+4703 obtained 
in the 1 .3 mm continuum and in the 12 CO(2- 1), 13 CO(2- 1), and 
C 18 0(2 - 1). Our observations at high angular resolution reveal 
bipolar outflow motion and gas accretion for IRAS 04579+4703. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations toward IRAS 04579+4703 were performed 
with the SMA on 2008 March 21, at 220 (lower sideband) 
and 230 GHz (upper sidebandfl- The two sidebands of the 
SMA covered frequency ranges of 219.4-221 .4 GHz and 229.4- 
231.4 GHz, respectively. The projected baselines ranged from 
7 k/l to 100 kA. The phase track center was R.A.(J2000.0)= 
05 h 01 m 39 s .92 and Dec.(J2000.0) = 47°07'21".10. The typical 
system temperature was 186 K. The spectral resolution is 0.406 
MHz, corresponding to a velocity resolution of 0.5 km s . The 
bright quasar 3C279 was used for bandpass calibration, while 
absolute flux density scales were determined from observations 
of Titan. QSO 0533+483 and QSO 0359+509 were observed 



1 The data are publicly available on the 
http://www.cfa.harvard.edu/rtdc/data/search.html 



website 
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Fig. 1. 1.3 mm continuum emission (black thicker contours) 
overlaid with the velocity integrated intensity map of 12 CO(2- 1) 
outflow (red and blue contours). The black contour levels are at 
3, 6, 9, 12, 15, 18 and 21 a- (1 cr is 0.002 Jy beam-'). The red 
and blue contour levels are 35, 50,..., 95 % of the peak value. 
The synthesized beam is 3".49 x 2".63, P.A.= -58.9° (shown in 
the lower right corner). The black plus symbol marks the peak 
position of the continuum source. The blue and red plus symbols 
represent the positions of the extracted spectra in Fig.4, respec- 
tively. H2O maser is shown with the filled triangle. The black and 
red open star symbols indicate the location of the infrared source 
IRAS 04579+4703 and a near-IR source detected by Varricatt et 
al. (2010), respectively. The blue and red dashed lines mark the 
direction of the NW and SE lobes of the H2 knots (Varricatt et 
al. 120 101) . respectively. 



al. I2008I I. Sanchez-Monge et al. (2008) argued that these radio 
emission could be from either radio jets or hypercompart Hn re- 
gions. The rising spectral index (Sanchez-Monge et al. 2008), 
association with an H2O maser, and compact source structure 
suggest a hypercompact Hn region in this region (Sewilo et al. 
2004). Therefore the 1.3 mm continuum emission may com- 
prise free-free emission from the hypercompart Hn region and 
the thermal emission from the warm dust from its surroundings 
(Keto et al. [19881 & Beltran et al. [20TTb . The emission from the 
free-free component of the continuum flux in our source can be 
estimated by S v oc y", where the free-free spectral index a is 1.1 
from 3.6 to 1.3 cm continuum emission. The contribution of the 
free-free emission to the 1.3 mm continuum is ~ 4.4 mJy (much 
lower than our total integrated flux 0.1 1 Jy at 1.3 mm), and then 
can safely be ignored. 

In order to determine the total gas mass of the continuum 
source, we make the optically thin approximation for the dust 
continuum emission. The total gas mass of the continuum source 
can be estimated by M core = S v D 2 /K v RB v (Td), where S ,, is the 
flux at the frequency v and D is the distance to the source. We as- 
sume k v — 0.05(v/230)^ m 2 kg -1 for the opacity (Testi & Sargent 
I1998K where we adopt/? = 2 (Drain & Lee 1 19841 ) and a dust-to- 
gas mass ratio R = 0.01 (Lis et al. I19911 >. B v (T d ) is the Planck 
function for a dust temperature Tj at frequency v. Using a dust 
temperature of ~ 30 K (Sanchez-Monge et al. 120081 ), we derive 
total gas mass = 13 M Q . Because the extended envelope can be 
filtered by the interferometers (missing short spacing data), the 
mass of the gas derived from our data is smaller than 23 M Q from 
the single dish observations (Sanchez-Monge et al. 2008). 



for antenna gain corrections. The calibration and imaging were 
performed in Miriad. A continuum image was constructed from 
the line-free channels. The spectral cubes were constructed using 
the continuum-subtracted spectral channels. Self-calibration was 
performed to the continuum data. The gain solutions from the 
continuum were applied to the line data. The synthesized beam 
size of the continuum was approximately 3".49 x 2".63 with a 
P.A. =-58° .9. Based on the phase monitoring of QSO 0533+483, 
the absolution position accuracy (pointing error) is estimated to 
be ~0.1" during the observations. 

3. RESULTS 

3.1. Continuum Emission at 1.3 mm 

In Figure 1, the 1.3 mm continuum image observed with 
the SMA shows a flattened source structure. Using a two- 
dimensional Gaussian fit to the continuum, we obtained a total 
integrated flux density of 0. 1 1 Jy, a deconvolved source size of 
4.7" x 2.9" (P.A. = 83.1°), and a peak position of R.A.(J2000) 
= 05 h 01 m 40. s 086, Dec.(J2000) = +47°07'20."90. The positional 
uncertainty due to noise can be estimated by A0 - 0.45 ^J™ , 
where 6*fwhm is the angular resolution and S/N is the signal-to- 
noise ratio (Rei d et al. 1 19881 Chen et al. 120061 Weintroub et al. 
2008; Qin et al. 120 101 ). The estimated positional uncertainty is 
approximately 0."14. If we also consider the uncertainty from 
the Gaussian fit to the position of the continuum peak and the 
pointing error, then AR.A. = +0.24" and ADec. = +0.21". The 
1.3 mm continuum emission is associated with an H2O maser 
(Migenes et al. 1999) as well as a near-IR source detected by 
Varricatt et al. (2010). In addition, the peak of the 1.3 mm con- 
tinuum agrees well with those of the 1.3 cm, 3.6 cm and 7 mm 
emission within the positional uncertainty (Sanchez-Monge et 



3.2. Molecular Line Emission 

Figure 2 shows the spectra of 12 CO(2 - 1), 13 CO(2 - 1) and 
C 18 0(2 - 1) emission at the peak position of the continuum im- 
age. The 12 CO(2 - 1) and 13 CO(2 - 1) lines present asymmetric 
profile with double peaks which can be caused by large optical 
depth (Andre et al. 1999 & Belloche et al . 120021 . The C ls O(2-l) 
line with a single emission peak can be considered as optically 
thin; Hence it can be used for determination of the systemic ve- 
locity. We measure a systemic velocity of 17.0 km s -1 from 

this line. The velocities of the blue and red emission peaks of 
12 CO(2 - 1) are -20 km s" 1 and -13 km s" 1 . For 13 CO(2 - 1), 
the velocities of the blue and red emission peaks are -19 km s -1 
and -16 km s . 

The 12 CO(2 - l) and 13 CO(2 - 1) spectra show "blue profile" 
signature (Wu et al. 2007), and the blue peaks of both lines are 
stronger than their red peaks with an absorption dip near the sys- 
temic velocity, which can be produced by infall motion or accre- 
tion. In addition, the I2 CO(2 - 1) spectrum show an unusual full 
width of ~30 km s -1 , from -38 to -8 km s -1 . Such a broad wing 
is a strong indication of outflow motion. 

To determine the velocity components and spatial extension 
of the outflow, we have made position-velocity (PV) diagram 
along the cuts at various position angles, and found that the red- 
shifted and blueshifted components are the strongest at a posi- 
tion angle of -40° (40° west of north) as shown in Figure 3. 
The PV diagram in Fig. 3 clearly shows bipolar components. 
The blueshifted and redshifted components have obvious veloc- 
ity gradients from -24.2 to -20 km s -1 and -13 to -8.5 km s -1 
respectively (see the blue and red vertical dashed lines in Fig. 3). 
The distributions of redshifted and blueshifted velocity compo- 
nents in Fig. 3 provide us further evidence for the bipolar outflow 
in this region (Yamashita et al. I19891 ). 
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Table 1. The physical parameters of the outflow. 

(xlO^Moyr- 1 ) 



Wing Velocity interval 
km s -1 



Deconvolved sizes 
arcsec x arcsec 



K 



T 
1 ex 

K 



N H2 

(xl0 21 cnr 2 ) 



M out 
(M ) 



td 
(xl0 4 yr) 



Blue -24.2 to -20.0 9.9x4.3(P.A.=-60.5°) 10.0 15.2 3.7 3.1 0.4 1.9 2.1 

Red -13.0 to -8.5 13.3x5.2(P.A.=-22.5°) 7.8 12.8 3.3 6.5 1.4 1.5 9.3 




Velocity (km s ) 

Fig. 2. Spectral profiles observed at the central position of IRAS 
04579+4703 in the optically thick 12 CO(2- 1), 13 CO(2- 1) lines, 
and the optically thin C ls O(2 - 1) line. The dotted line in the 
spectra marks the cloud systemic velocity. 



Using the velocity ranges determined from the PV diagram, 
we made the velocity integrated intensity maps superimposed on 
the 1.3 mm continuum as shown in Fig. 1. The blueshifted and 
redshifted components are presented as blue and red contours. 
The map clearly shows northwest-southeast bipolar components 
centered at the peak position of the 1.3 mm continuum. The 
blueshifted and redshifted lobes are not in a straight line. The 
outflow is associated with infrared source IRAS 04579+4703 
and a near-IR source. 



4. DISCUSSION 

4.1. Outflow Kinematics and Driving Source 

The velocity integrated intensity map, PV diagram and the broad 
wing (full width = 30 km s -1 ) of the CO 12 CO(2 - 1) line 
clearly show a bipolar outflow in IRAS 04579+4703 region. 
The blueshifted and redshifted lobes do not lie in a straight line, 
and the extension is similar to that of H2 knots (Varricatt et al. 
1201 Oi l, indicating that the H2 and CO most likely represent the 
same outflow, and that the 12 CO(2 - 1) outflow gas probably 




Fig. 3. The position- velocity plot of the 12 CO(2 - 1) emission 
at a position angle of -40°. The contour levels are 1.4, 1.8, 2.3, 
2.7, 3.2, 3.6, 4.1 and 4.5 Jy beanr 1 . 1 cr is 0.7 Jy beam" 1 . The 
black vertical dashed line marks the cloud systemic velocity. The 
blue and red vertical dashed lines mark the velocity ranges of the 
blueshifted and redshifted emission, respectively. 



is entrained by the H2 jet-like outflow. The outflow is associ- 
ated with the infrared source IRAS 04579+4703 and a near-IR 
source. Varricatt et al. (2010) suggested that the near-IR source 
may be the near-IR counterpart of IRAS 04579+4703 with the 
strong infrared excess. The position of the near-IR source agrees 
well with the peaks of the 1.3 mm, 7 mm, 1.3 cm and 3.6 cm 
emission within the positional uncertainty. H2 and Bry emission 
are very close to the near-IR source (Ishii et al. 2001 & Varricatt 
et al. 120101 1. Together, we conclude that the near-IR source may 
be an embedded massive YSO, which is likely to drive the bipo- 
lar outflow. 

Under conditions of local thermodynamical equilibrium 
(LTE), 12 CO(2-l) being optically thick, and the completely fill- 
ing beam with the source (the filling factor / = 1), the column 
density of the outflow is estimated by using the equation (eq. 1) 
following Garden et al. Q991I I with a relation JV H , ~ l0 4 Nn co ( 
Dickman fT978l : 



N Hl = 1.08x10' 



(r ex + 0.92) 
exp(-16.62/T ex ) 



/ 



T mh x 



rdv 



[l-exp(-r)] 



cnT 2 ,(l) 



where T ex is the excitation temperature, T m \, is the 
corrected main beam temperature, and r is the optical 
depth for the 12 CO(2-l) line. We assume the solar abun- 
dan ce rat io [ 12 CO]/[C 18 O]=t( 12 CO)/t(C 18 O)=490 (Garden et 
al. [1991). T ex is estimated following the equation r ex = 
ll.l/ln[l + 1/(7^/11.1 + 0.02)] (Garden et al. [19911 The red 
and blue lobes of the outflow peak near the continuum, then we 
assume that the 12 CO and C ls O trace same cloud components in 
the peak positions of the two lobes for the calculation, the as- 
sumption is reasonable since both I2 CO(2 - 1) and C ls O(2 - 1) 
are low energy transition with similar upper level energy. Under 
assumption of C 18 being optically thin, we calculate t accord- 
ing to r = -4901n[l - (e lu/7 « - 1)7^/11.1] (Garden et al. 
1991). Under the Rayleigh-Jeans approximation, 1 Jy beam 1 
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Fig. 4. C 18 0(2- 1) spectra. In a and b panels, the molecular spec- 
tra are extracted from the positions of the blue and red plus sym- 
bols marked in Fig. 1, respectively, which are separated by 4". 
The black vertical dashed line marks the cloud systemic velocity. 
The blue and red vertical dashed lines mark the peak velocities, 
respectively. 

in our SMA observations corresponds to a brightness tempera- 
ture of 2.52 K. Furthermore, the mass of the outflow is given 
by M out = mNn 2 S/(2.0 xlO 33 ) M , where m is 1.36 times the 
H2 mass (Garden et al. 1991) and S is the size of each lobe. S 
is obtained from the deconvolved lobe size (In Table 1). A dy- 
namic time scale can be determined by td - r/v, where v is the 
maximum flow velocity relative to the cloud systemic velocity , 
and r is the length of the begin-to-end flow extension for each 
lobe. The lengths of the blueshifted and redshifted lobe are 0.14 
pc and 0. 1 3 pc obtained from Fig. 1 , respectively, which is only a 
lower limit due to projection on the sky plane. The mass entrain- 
ment rate of the outflow is calculated using M oul - M out /(fd). 
The physical parameters and the calculated results of the out- 
flow are summarized in Table 1 . The parameter difference of the 
two lobes reflects that the ambient material is not uniform and 
then the entrained masses and dynamics are different. The out- 
flow has the total mass of 1.8 M G and the mass entrainment rate 
of 1 . 1 x 1CT 4 M Q yr~ l . The average dynamical timescale is about 
1.7xl0 4 yr. 



4.2. Accretion 

Our interferometric continuum image toward IRAS 04579+4703 
reveals a flattened structure, perpendicular to the direction of the 
outflow detected in 12 CO(2 - 1) and H2, which is consistent with 
accretion disk scenario observed in the other high mass star for- 
mation regions (Cesaroni et al. 1999; Shepherd et al. 2001; Patel 
et al. 2005). The C 18 0(2 - 1) spectra from the two positions (the 
red and blue plus symbols in Fig. 1) perpendicular to the direc- 
tion of the outflow are shown in Figure 4. The velocity gradient 
is seen from the C 18 0(2 - 1) spectra, suggesting that the disk 
is rotationally and gravitationally bound by the central YSO. 
Assuming the rotation is Keplerican, the dynamical mass (bind- 
ing mass) is estimated by M - 6v 2 r/Gsin 2 (i), where r is the half 
spatial separation of the emission, 6v is the difference between 
the peak velocity of the emission and the systemic velocity, and 
i is the unknown inclination angle between the disk plane and 
the plane of the sky (i = 90° for an edge-on system). Adopting 
the values r ~2" and 5v ~ 1 km s shown in Figures 1 and 4, 
the obtained dynamical mass is ~6/[sin 2 (O]M . The dynamical 
mass is estimated to be 12 M if i = 45°, which is a lower limit 
since that we estimated the binding mass from the kinematics 
inside 2" due to our low spatial resolution observations. 

The spectroscopic observations of H2 and Bry (Ishii et al. 
120011 ) also suggested the presence of an accretion disk, and 
this source is probably accreting. If we assume that the cen- 
tral star is a massive protostar and not in the phase of zero- 



age main-sequence (ZAMS) star, then the bolometric luminos- 
ity (3.91xlO 3 L ) is equal to accretion luminosity. According 
to Molinari et al. dl998bl l, the accretion rate is given by M = 
6.22 x 10- 9 (L/L o ) L70 (M/M G r L24 . Assuming the dynamical 
mass of 12 M is equal to the protostar mass, the expected ac- 
cretion rate is 3.7xl0~ 4 Af yr _1 . 

The "blue profile" of the 12 CO(2 - 1) line provide a further 
evidence for infall motion or accretion as discussed in section 
3.2. The mass accretion rate is estimated by M - AjiR 2 Vmn 
(Myers et al. 1996), where m is 1.36 times the H2 mass. If cloud 
core is approximately spherical in shape, the mean number den- 
sity is n- 1.62 x 10~ 19 /VH,/i-, where L is the cloud core diameter 
in parsecs, which is 0.056 pc (the major axis of the deconvolved 
size of the continuum). N(H2) is the column density of molec- 
ular hydrogen obtained from the equation (eq. 3) of Xu et al. 
d20T0b in C ls O line. The velocity difference (V) of 2.0 km s _1 
between the systemic velocity (-17.0 km s -1 ) and the velocity of 
the redshifted absorbing dip (-15.0km s -1 ) and a source size (R) 
of 0.028 pc, then we obtain a mass accretion rate of ~ 2.3xl0 -4 
M yr~' which is roughly consistent with the value derived from 
bolometric luminosity, further suggesting the presence of an ac- 
cretion disk in this region, but the observations at the current 
spatial resolution are not adequate to distinguish the disk accre- 
tion and the observations of more molecular tracer with higher 
spatial resolutions are needed to resolve detailed kinematics in 
this region. 
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